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We present a theory of spin waves in a non-condensed gas of spin-1 bosons: providing both 
analytic calculations of the linear theory, and full numerical simulations of the nonlinear response. 
We highlight the role of spin-dependent contact interactions in the dynamics of a thermal gas. 
Although these interactions are small compared to the thermal energy, they set the scale for low 
energy long wavelength spin waves. In particular, we find that the polar state of 87 Rb is unstable 
to collisional mixing of magnetic sublevels even in the normal state. We augment our analytic 
calculations by providing full numerical simulations of a trapped gas, explicitly demonstrating this 
instability. Further we show that for strong enough anti-ferromagnetic interactions, the polar gas is 
unstable. Finally we explore coherent population dynamics in a collisionless transversely polarized 
gas. 



I. INTRODUCTION 

To date, all spin wave experiments in non-condensed 
systems have been on spin-^ or pseudo-spin-^ systems. 
Here we ask the following questions: what is the nature 
of spin waves in a thermal spin-1 gas? How do these differ 
from the well understood spin-^ case? These questions 
are of experimental interest, with several groups possess- 
ingthe technology to study them using ultra-cold gases 
III- 

In the context of cold gases, spin waves were first dis- 
cussed by Bashkin and independently by Lhuillicr 
and Laloe @. The key finding was that spin exchange 
collisions can give rise to weakly damped spin waves even 
in a non-degenerate gas. Ultra cold gases have provided 
an exciting setting for observing these spin phenomena. 
In particular, experiments on pseudo spin-i Bose and 
Fermi systems by the JILA @ and Duke Q groups have 
observed coherent collective oscillations in an otherwise 
classical gas. 

One expects the physics of a spin-1 gas to be far 
richer than a pseudo-spin- \ system. This is amply 
demonstrated by experiments on condensed spin-1 gases 
B H, EH EJ- One of the most dramatic observations 
was that of a dynamical instability in Bose condensed 
87 Rb, studied by Sadler et al. Beginning with a gas 
pumped into the ttlf = state, they observed the spon- 
taneous formation of transverse ferromagnetic domains. 
Will a similar instability be observed in the normal state? 
We find that there is an exponentially growing mode in 
the unmagnetized gas, even in the normal state, with a 
wavelength comparable to typical cloud sizes. 

From a theoretical perspective, the source of novel 
physics in a dilute spin-1 Bose gas (such as 87 Rb, 23 Na) 
is the structure of the interactions, described by two cou- 
pling constants: cq and C2, representing spin-independent 
and spin-dependent contact interactions. The intcrac- 



* Electronic address: ssn8@cornel l.edul 
^Electronic address: cm256@corncll.edu 



tion Hamiltonian density takes the form 77; nt = cqti 2 /2 + 
C2(S) 2 /2, where n and S are the local density and spin 
density respectively [l3[. The coefficient C2 has no cor- 
responding analog in the spin-^ case. This interaction 
gives rise to spin mixing collisions, where two atoms in 
the mjf = hypcrfine sub-level can combine to form 
atoms in the mp = ±1 states 1J]. Another important 
consideration is the quadratic Zeeman effect, which arises 
from the hyperfine interaction and the difference in the 
coupling between the electronic and nuclear spins. The 
linear Zeeman effect can be neglected in the Hamiltonian 
as the total spin of the atoms is a conserved quantity. 



We begin our analysis by setting up the problem, and 
reviewing spin waves in a spin-i gas. This allows us 
to highlight the differences with the spin-1 case. Next, 
we turn to the spin-1 gas. Starting from a microscopic 
Hamiltonian, we obtain a linearized Boltzmann equation 
about the ferromagnetic (mj = 1) and polar (rriF = 0) 
states, and calculate the spin wave dispersion in each 
case. We find that a polar gas with ferromagnetic in- 
teractions is dynamically unstable towards spin-mixing 
collisions for small enough Zeeman fields, analogous to 
the condensed case. By explicit calculation of the dis- 
persion relation, we show that for strong enough anti- 
ferromagnetic interactions, an instability occurs in the 
polar state [T^. 



Following these analytic calculations, we perform nu- 
merical simulations of a trapped gas using an effective 
ID Boltzmann equation. We explore the evolution of 
transversely polarized spins, and investigate dynamical 
instabilities. 



Our work complements prior work on the kinetics of a 
normal spin-1 Bose gas by Endo and Nikuni [2]. While 
their focus is on the effect of collisions in the damping 
of collective modes of a trapped gas, with particular em- 
phasis on dipole modes, we focus here on the collisionless, 
Knudscn regime. 
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II. BASIC SETUP 



A. Kinetic Equations 



As described by Ho [HI and Ohmi and Machida [13], 
the second-quantized Hamiltonian for a spin-1 gas, ex- 



H= / dr^ a 



~2^ + U ^ 



qS z - S z ip, 



pressed in a frame where each spin component is rotating 
at its Larmor frequency is 



r 



(i) 



where a = ( — 1,0,1), is the quantum number for the z- 
component of the spin, Vv( r ) is the field annihilation 
operator obeying bosonic commutation relations, and 
U(r,t) is the trapping potential. Here S denotes the di- 
mensionless vector spin operator. Throughout, we set 
h = 1. Boldface is used to denote matrices, and arrows 
denote vectors. 

The interaction strengths, expressed in terms of the 
scattering lengths in the spin-2 and channel (ct2, do) 
are Co = 47r(ao + 2a2)/3m and C2 = 47r(a2 — ao)/3m. A 
negative C2 favors a ferromagnetic state with (S) = 1, 
while for positive C2, the equilibrium state is one with 
(S) = 0, where all the atoms are in the tuf = state, or 
simply an incoherent mixture. 

Additionally, one considers the quadratic Zeeman ef- 
fect q oc B 2 , which favors a state with tuf = ((S- S) = 
0). Thus for a gas with negative C2 such as 87 Rb, the 
spin dependent contact interaction competes with the 
quadratic Zeeman field, giving rise to interesting dynam- 
ics [3. 

In the condensed gas, dipolar interactions may also 
be important [H, Q . At the lower densities found in a 
normal gas these interactions, which fall as 1/r 3 , may be 
neglected. 

Following standard arguments [l8[ we obtain the 
equations of motion for the Wigner density ma- 
trix F a b(p,R,t), whose elements are f a b(p,R,t) = 
J dre-^ f {rpl(R - ^,t)tp b (R + f,t)>. The diagonal com- 
ponents of the spin density matrix, when integrated in 
momentum, give the densities of each of the spin species. 
The off diagonal terms, often referred to as coherences, 
are responsible for spin dynamics. Here p represents mo- 
mentum, R and r denote the center of mass and relative 
coordinates. The Wigner function is the quantum analog 
of the classical distribution function. By taking moments 
of the Wigner function, we obtain physical observables 
such as the density n(R,t) — J j0yz F(p,R,t) and spin 

current t) = J p F(p, R 7 1). 

The equation of motion takes the standard form 0]: 



dt 



V,F 



{VrV,V p f}+I c (2) 



where I c is the collision integral and V is the interac- 
tion potential. This form of the Boltzmann equation (|2|) 
is completely general, and holds for any non-condensed 
spinor gas. The role of spin enters in determining the 
dimension of the density matrix, and the exact form of 
the interaction potential V. 

The first two terms on the right hand side of the Boltz- 
mann equation @ arise from forward and backward scat- 
tering collisions. While the former type merely alter the 
mean field seen by the atoms, backward scattering colli- 
sions allow the colliding atoms to exchange momentum. 

The last term in the Boltzmann equation is the colli- 
sion integral, responsible for energy relaxation. While a 
detailed derivation of the collision integral is non trivial 
fl6| , the qualitative properties are well described within a 
simple relaxation time approximation I c = — (/ — /o)/t. 
The relaxation time (r) is proportional to the elastic scat- 
tering rate rj 1 ~ 87raQU T n, where v T = \J 2k ^ is the 

thermal velocity and n is the density. 

The precise expression for V determined by including 
the interactions within the Hartree-Fock description is 

V=(U + c Tr(n)) 1 + c n + c 2 Sn- S + c 2 M- S (3) 

where 1 is the identity matrix, Tr denotes the trace op- 
eration, and M = Tr(S n) is the magnetization. Our 
form for the interaction energy is equivalent to that of 
Endo and Nikuni [ijj]. One may explicitly check that 
this interaction potential is rotationally invariant in spin 
space. 

In experiments the external trapping potential U is 
often spin-dependent, an effect which is readily incorpo- 
rated. The spin independent interaction gives rise to self 
and cross interaction terms. The latter give rise to co- 
herences and are encoded in the second term in In 
addition, the mean-field potential alters the trapping po- 
tential seen by all the atoms by an amount contot, where 
the total density n to t = Tr(n). 

The contribution of the spin-dependent interaction is 
more subtle and can be understood as follows: The first 
term accounts for spin dynamics such as spin-relaxation 
collisions. For ferromagnetic interactions (02 < 0), the 
last term increases the density of regions where the 
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atomic spins arc aligned with respect to one another 
(|M| = 1). For a fully polarized gas in the m F = +1 sub- 
level, this is C2S2, while for an unmagnetized {mp = 0) 
gas it is zero. 

Finally, we note that while the self- and cross- interac- 
tions between the three sublevels produce diagonal and 
off-diagonal contributions to the interaction potential, 
the spin-relaxation or population exchange interactions 
only gives rise to coherences, and are absent for a single 
component gas. 



B. Qualitative Features 

Here we elaborate on our argument for why one ex- 
pects richer spin physics in a spin-1 system versus a 
spin-i one. We assume a uniform, collisionless gas with 
V RCon/mvT <C cqji, such that all the physics is governed 
by the commutator in ([2]). 

The pseudo-spin ^ case may be understood starting 
from the fact that that a 2 x 2 matrix (A) can be de- 
composed into A = A 1 + A^a^, where A = Tr(A), and 
a arc the Pauli matrices. Expressing the density matrix 
and interaction potential in this way, the equations of 
motion are §- t F Q + |- V fl F = and f + 1- = 
~FpV v t pini , where e pl/M is the completely antisymmetric 
unit tensor. The second equation, which is responsible 
for much of the spin wave physics in a spin- i system sim- 
ply says that the interactions act as an effective magnetic 
field about which the spins precess. 

To anal yze the spin-1 case, we follow Ohmi and 
Machida [17| and use a Cartesian basis (ip = 
{ip x , ipy, ipz})- This representation is related to the spher- 
ical ({1,0,-1}) basis as follows [H]: ip x = "^(V'l — 
ip-i),ip x = -^=(ip 1 + tp_x) and ip z = ip . 

In the Cartesian basis the irreducible decomposition 
of a spin-1 system is A = AqI + ie^cA^ + A^ s \ where 
the scalar Aq = , A^ = e a t> c A^ is a completely 
antisymmetric matrix proportional to the vector spin (S) 
order, and A' s ' is a symmetric traceless tensor which is 
related to the spin fluctuations (S a Sj, + S;,S a ), and is a 
nematic degree of freedom. 

Writing the density matrix (F) and interaction (V) in 
terms of their respective irreducible decompositions we 
obtain three equations of motion: Jj-Fo + ™' ^ R-^o = 0, 
JLpto + £..y R F^) = l [V (Q) ,F (a) ] -M[V™ s) ,F (s) ], and 
9_F^ + l^jjfW = i[V (a) ,F (s) ] + i[V^,F (a) ]. Here 
the second equation describes the evolution of (S) and the 
spin currents, while the last equation describes dynamics 
of the nematicity (A^). 

The new feature is that the spin and nematic degrees of 
freedom are coupled by the interaction matrix. Using the 
Levi-Civita identities, it is easy to see that i[V^ a ' ,F^ a ') 
is the spin-1 analog of the corresponding spin-i term 
which acts as an effective magnetic field for the spins. 
But the last term i[V^ s \ F^] shows that fluctuations of 



the nematicity can change the spin dynamics. 



C. Spin waves in a two-component gas 

As a starting point for understanding the spin-1 gas, we 
review the spin-i case. As explained in [20| . spin dynam- 
ics in the cold collision regime are governed by momen- 
tum exchange collisions with a characteristic timescale 
T ar 1 ~ 47r °^ 2 " , where an, is the two-body scattering 
length. When ksT <C m ^ , or alternatively, the thermal 
deBroglie wavelength is large compared to the scattering 

length A T = \J ,nk B T — ai2 > Tex ^ Tei ' and for timcs 
shorter than r e ;, particles exchange momentum several 
times without significantly altering the energy distribu- 
tion. Therefore excitations with wavelength longer than 
VTT ex are a collective effect. 

For typical densities (~ 10 14 cm -3 ) and scattering 
lengths, n 1 / 3 a ~ 0.01 the condensation temperature (T c ) 
for the interacting system may be approximated to that 
of an ideal Bose gas [2l|. Thus in order to see collective 
spin phenomena, we require T c ~ (lo)N^ < ksT < , 

where (lu) = ( s L0 x L0 y L0 7 )) 1 l' i is the average trap frequency 
[22| . Although this is typically a wide temperature range 
10~ 7 < T < 10~ 3 K, diffusive relaxation damps out spin 
waves at higher temperatures. 

A two-component gas has a longitudinal and a trans- 
verse spin mode. The longitudinal mode has a linear 
dispersion, but is strongly Landau damped By con- 
trast, the transverse mode is weakly damped and prop- 

agates with dispersion lo = — ^ Tex f 1 — i ^ J , where 

~ < 1 [23j ]. These weakly damped spin waves 
have been observed in the context of spin polarized Hy- 
drogen by Johnson et al. @, 0, HI] • 



III. SPIN WAVES IN THE SPIN-1 GAS 

In this section we linearize Eq.([2]) about stationary 
states. In Sec. IV wc consider more general dynamics. 



A. Excitations about the mp — 1 state 

We begin by linearizing about the niF = 1 state. 
We start with a homogenous gas of particles with a 
Maxwcllian velocity distribution and initial density tiq. 
The self-interaction between two particles is proportional 
to 0,2, and so we define fij„t = (cq + 02)710. We begin by 
considering the collisionless limit, 1 <C Oj„ t TD. 

Linearizing about this state, and dropping the collision 
term, the Boltzmann equation can be written in Fourier 
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FIG. 1: (Color Online) (a) Dispersion relation ui 2 (k) for var- 
ious values of q using 87 Rb parameters (02 < 0) at T = luK 
and no = 10 14 cm -3 : q=0 (solid, black), q = 2|c2|no (blue, 
dashed), q = 4|c2|no (green, dotted), q = 4.5|c2|no (red, thin), 
(b) Gapless, linear dispersion for small k at q — 4\c2\uq. 
The solid black curve is the small k expansion J7J. (c) 
q = 2jc2|no dispersion in detail. The horizontal black line 
is uj 2 — g(q — 4|c2|no), which is the gap predicted by the small 
k expansion (J7J. (d) q — dispersion in detail. Horizon- 
tal dashed line at ui 2 = —Ac^n 2 indicates the most unstable 

mode. Vertical dotted line is at kAr = V 2 ^ 2 ^ " ".. 



space: 



TO 



Vn,(5F 



^V,F j 
fc <5V,V P F } 



(4) 



where Fo is the initial distribution and Vq is the ini- 
tial interaction potential. The quantities <5F and SV 
are Fourier transforms of the change in the distribution 
and interaction, defined in the usual way 5f(k,p,u>) = 
J drdte l ( k ' r ~ UJt ' > 8f(p,r,t). Inserting the form for the in- 
teraction potential ([3]) , we obtain the equations of motion 
for the density and spin fluctuations. 

We find that four long wavelength (|fc|t^r <C Hint), l° w 
frequency (u <C Hint), spin modes propagate in this sys- 
tem. These modes have a quadratic dispersion, charac- 
teristic of ferromagnetic systems, but have a k = off- 
set. Two of the modes are transverse spin waves, and the 
other two are quadrupolar modes where the longitudinal 
magnetization fluctuates along with the transverse spin 
fluctuations ({S»S„) - (S^S^)). 

The spin modes are given by the equation 1 = — (co + 
„ \ r dp /o+'jfc-Vp/o 

C 2 )J Jz^S 



uj+n z „ t +q-^ 

have dispersion uj(k) = 



25j. The transverse spin modes 



■-!<>,, 



— q. The fact that this has a 



negative energy at k = signifies that in the presence of 
a quadratic Zeeman shift, the initial state is not the ther- 
modynamic ground state. One can decrease the energy 
by rotating the spins to be transverse to the magnetic 
field. However the frequency stays real, meaning that 



the initial state is mctastablc, and transverse spin will 
not be spontaneously generated. 

Meanwhile, the quadrupolar modes have dispersion 

uj(k) = —4c2Uq + 2n T t . In a gas with positive C2, the 
frequency vanishes at a finite wave-vector, indicating a 
thermodynamic instability in the anti-ferromagnetic gas. 
Conversely for negative c 2 the spin-dependent contact 
interaction leads to a gap in the quadrupolar mode spec- 
trum. 

As pointed out by Leggett [2(|, one can also derive 
these equations by writing the macroscopic equations for 
the density and spin current: 



d t n 



V R -j 



V,n 



(5) 



dtj 



V R Q 



VJ 



-5{*»m4 



where the energy density Q = J dp (p-p) F. In order 
to obtain a closed set of equations, we approximate Q w 
J dp | (p-p) F = ^n, where pt = ttivt- In a collisionless 
gas, sufficiently close to equilibrium, these approaches are 
equivalent. 

An advantage of the Leggett approach is that it pro- 
vides access to the dispersion at large k more readily 
than the former approach. Solving (J5), we get the stan- 
dard relations u(k) = q - h+fn^J^s) (^mtTo - i), and 

Lu(k) = -4c 2 n + (i+(o^p 2 ) 2 ) (HintTD " 0- Setting td to 
infinity, we once again recover the relationships derived 
above. 



B. Excitations about the mp — state 

Next we consider the case of a gas in the mp = 
state. As before, low energy dynamics are driven by a 
combination of the quadratic Zeeman energy, and the 
spin-dependent contact interaction. As we now show, in 
the case where c 2 is negative, these energies compete, 
driving a dynamic instability. Linearizing equation (j4|) 
yields the relation: 



1 



4C2X1X2 



(1 + (c + c 2 )x 2 )(l - (co + c 2 )xi) 



(6) 



where the response 

r dp (fo + ik-V p J ) 



and X2 



functions are xi 

' dp (/o-^fc-Vp/o) 



and Orf = (co — C2)?io — Q- For \k\vx "C cqUo and 
lu <C cqUq. The resulting dispersion is: 



u 2 (k) = q(q + 4c 2 7i ) + (q + 2c 2 n ) 



k 2 v 2 T 
c n 



(7) 



The spin and quadrupole modes in the thf = state 
are degenerate. For C2 > 0, or for large enough q, one ob- 
tains a gapped quadratic mode. For c 2 < however, the 
dispersion goes from quadratic to linear at q = 4|c2|no. 
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FIG. 2: Microwave-field induced instability in a polar gas for 
q < 0: The modes of a polar gas with ferromagnetic interac- 
tions (c2 < 0) are plotted for q = Acin (top) and q = 2c2n 
(bottom) (hence q < 0). The gaps are given by q(q + 2c2«.o) 
in each case. Note that for negative q, a finite wave- vector in- 
stability sets in. Increasing \q\ pushes the instability to larger 
wave- vectors. The dashed line is the most unstable frequency 
for q = (see Fig. [IJd)) given by uj 2 = -(2c 2 n ) 2 . The 
timescale for the onset of the instability is independent of q. 
We set T = 1/xK and n = 10 14 cm -3 . 



2 i 2 

((w - n a )(w + n d ) - (( w + n.)(w - n, 



where f2 s = g + 2c2no- Expanding ((Sj) to lowest order in 
k returns Eq.0. 

The roots of (|5J) for negative C2 are plotted in Fig.QJa) 
for various values of the quadratic Zeeman energy q. 
The temperature is set to 1/tK, density to 10 14 cm" 3 , 
and we use the scattering lengths for 87 Rb, ao(a2) = 
5.39(5.31)nm. For the homogenous gas in consideration 
here, this is well above the temperature for Bose-Einstein 
condensation. As shown in Fig. |Tfb) , at q = 0, u> 2 (k) is 
negative for small k, indicative of an instability. 

The timescale for the onset of this instability (U ns ) is 
determined by the most unstable mode, which from © 
occurs at u? = -(2|c 2 |n ) 2 , giving t ms ~ 2 \cl\n - The 

corresponding wave- vector is kinsA-T = ^ 2 |^y°" -. Solv- 
ing for where the mode frequency vanishes, we find that 
all wave-vectors k smaller than a critical wave-vector 
k c = \ 2ki ns lead to unstable modes. For achievable ex- 
perimental densities (~ 10 14 cm~ 3 ) and temperatures (~ 
1/xK), this wave-length of the instability A c ~ lOOA^ = 
10/tm, which is comparable to typical cloud sizes. The 
time for the onset of this instability ~ 0.25s. Although 
this timescale is easily achievable in experiment, small 
collisional timescalcs (~ 50ms) at these temperatures 
may make the instability undetectable in experiment. 

Microwave induced instability in a ferromagnetic polar 
gas:— [32j Gerbier et al. (27j have shown that a weak, 
microwave driving field, off resonant with the F = 1 — > 
F = 2 transition, may be used to tune the energy levels 
of the F = 1 hypcrfine sublevcls. In this sense, it plays 
the same role as the quadratic Zccman shift. By changing 



Upon lowering q further, the rap = state undergoes 
a dynamic instability where the excitation frequency is 
complex. At q = 2|c2|no, the leading order term in the 
small k expansion vanishes. 



The structure of the modes should be contrasted with 
those in the mj? = 1 state. Even though either case 
may not be a thermodynamically stable starting point 
for anti-/ferromagnetic interactions (respectively), only 
in the antifcrromagnctic case do small fluctuations grow 
exponentially. 



Taking the tjj — > oo limit and solving (O, the disper- 
sion relation is determined by the equation: 



)-^-) =-4(c 2 n ) 2 (w-fi d )(u; + fi d ) (8) 
I 

the detuning, one can change the effective q from positive 
(E mF = < E mF=±1 ) to negative {E mF=0 > E mF=±1 ), 
where E denotes the energy of the hyperfine sublevels. 

It is therefore of experimental and theoretical interest 
to consider what happens to the rap = 1,0 states at neg- 
ative q. For the rnp = 1 state, the spin mode becomes 
gapped by the quadratic Zeeman energy, and the ther- 
modynamic instability disappears. This is to be expected 
as a negative q favors a ground state with (S- S) = 1. 

The polar state however is more interesting. From l[7|). 
it is clear that for ferromagnetic interactions (c2 < 0), the 
first term on the right hand side is always positive, and 
spectrum becomes gapped. The second term however is 
negative, signaling a finite wave-vector instability. Ana- 
lyzing © (sec FiglJ), we find that increasing \q\ pushes 
the instability to larger wave-vectors, but frequency of 
the most unstable mode remains largely unchanged. The 
experimental consequence is that q < 0, one should ob- 
serve a reduction in the size of ferromagnetic domains, 
while leaving the timescale for the onset of the instabil- 
ity unaffected. 

Large anti-ferromagnetic interactions: The modes of 
the polar gas, when c 2 > can be obtained from (jSJ. 
For small C2/C0, the polar state is stable toward spin and 
nematic fluctuations. However, when c-2 and cq are com- 
mensurate, an instability can set in. This latter limit, 
which has been largely unexplored, is primarily of aca- 
demic interest as all current spinor gas experiments have 
I C2 1 *C Co. For simplicity we restrict our analysis to the 
case of no magnetic field (q = 0). 
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FIG. 3: (Color Online) Dispersion relations (uj(k)) for C2 > and ~ Co. (a) Gapped (top) and ungapped (bottom) quadratic 
dispersion relation for C2 = 0. (b) The two solid curves in (a) are now shown for C2 = 0.5co At some wave-vector the two 
modes become one single mode. The real (black) and imaginary (orange) components are indicated separately. The gap is 
equal to fid = Co — C2. (c) Linear dispersion uj{k) = -j^vrk, when C2/C0 = 1. (d) Two modes for C2 = 1.5co: Gapped quadratic 
real mode; an unstable mode: real (black) and imaginary (orange) parts shown. Notice that the instability only occurs if 
< k < kmaxipv/co). The parameters for these plots are T = 1/iK, n ~ 10 14 cm -3 . For definiteness, we picked the atoms to 
have the mass of 87 Rb and the corresponding Co value. We caution however that this choice is artificial, as | C2 1 <C Co in 87 Rb. 
All these plots were made in zero field (q — 0). 



As shown in Fig. EJa), for C2 = we find a gapped 
and ungapped quadratic mode with dispersions ui{k) = 

J \{Q 2 d ± \/^ d ~ ( v Tk) 4 ), where fid reduces to cono in 
this limit. As ci is increased, for some k c , the gapped 
an ungapped modes meet, and for k > k c , imaginary fre- 
quencies appear (Fig. E^b)). While the exact expression 
for k c is complicated, and depends on T and density, 
we estimate the appearance of imaginary wave- vectors 
when fi s ?=a Q d i.e. C2 ~ co/3. The case of C2 = cq 
is special, and there is a single linear mode with fre- 
quency uj(k) = -^vxk (Fig. (3Jc)). However as C2 > c 

(Fig- Eld)), imaginary frequencies appear once again for 
a bounded set of k values: < k < k max (c2/co). 



C. Summary of Results 

In Fig. 0] we summarize the results of our stability 
analysis in terms of the interaction parameters cq and C2 
for zero q. The following general features are clear: 

1. For Co > and C2 < 0, the ground state is a ferro- 
magnet. 

2. For Co > and C2 > 0, but C2 < cq/3, the ground 
state is polar. 

3. In region III the scattering length in the spin-2 
channel (02) is negative, and the gas becomes me- 
chanically unstable. 

4. In region IV, the polar state is metastable, in the 
presence of finite populations of ±1 components, 
the system will phase separate [lj| . The ferromag- 
net is thermodynamically unstable here. 

5. The regime of current day experiments on spin-1 
gases is that of \c2\ <C cq, indicated by the open 
and filled circles in the figure. For 87 Rb, at 
10 14 cm" 3 , c n ~ 50nK and \c 2 \n - 0.25nK (filled). 
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FIG. 4: Stability phase portrait in terms of the interaction 
strengths con and C2n at zero Zeeman field. The interaction 
strengths are measured in nK to facilitate comparison with 
typical experimental values. The four distinct regions are de- 
marcated by solid black lines. I: For C2 < 0, but C2 > -co, 
the ferromagnetic state is the thermodynamic ground state. 
II: For C2 > 0, C2 < Co/3, the polar state is the ground state. 
Ill: Here a^{= co + C2) < 0, and the gas is dynamically un- 
stable towards collapse. IV: When C2 > Co/3, the polar state 
becomes metastable. In the presence of finite populations of 
±1, phase separation is expected to occur. In this region, the 
ferromagnetic state is thermodynamically unstable. Finally, 
the open and filled circles are typical interaction scales in 23 Na 
and 87 Rb for n ~ 10 14 cm -3 . 



For 23 Na, at the same density cqu ~ 85nK and 
C2n ~ 4.25nK (open). 



IV. NUMERICAL INVESTIGATION 

In an experiment, the atoms are typically confined to a 
harmonic trap. Calculating the modes in a trap is much 
more difficult than in free space. Example calculations 
[l6| do not always yield simple physical pictures. More- 
over experiments are typically performed in a limit where 
linear response is not applicable 0]. 

For these reasons, we perform a numerical study of the 
spin dynamics for a trapped gas. Furthermore while the 
calculations in part I pertained to the collisionlcss regime, 



using a simple relaxation time approximation, we model 
collisions in the gas. 

We address two questions here: 

1 . What happens if the initial state is not a stationary 
state with respect to the external magnetic field ? 

2. Can the instability in the polar state be observed 
experimentally? 

To address question 1, we consider a ferromagnet in 
the x direction, in the presence of a field along z. We 
find that of coherent population dynamics in a collision- 
less gas. By controlling the rate of decay of spin current, 
one may be able to experimentally observe coherent os- 
cillations even in a classical spinor gas. 

With regards to 2, taking into account realistic experi- 
mental parameters, the polar state instability should not 
be visible in 87 Rb. The collision time is set by Co, while 
the time for spin dynamics by C2 . Owing to the small ra- 
tio of | C2 1 / Co ~ 0.005, collisions lead to relaxation before 
any coherences can develop. 



A. Numerical setup 

We consider parameters such that the all the mo- 
tion takes place in one dimension producing an effective 
ID Boltzmann equation. Assuming that the distribu- 
tion function in the transverse directions is frozen into a 
Boltzmann form, we can integrate them out [30l |. yield- 
ing renormalized scattering lengths (00/02). In what fol- 
lows, we normalize all lengths by the oscillator length 
y/l/(mu} z ) and momenta by ^/rnksT. Using a phase- 
space conserving split-step approach (3lj | , we numerically 
integrate @. We have verified that the tempero-spatial 
grid is fine enough that our results are independent of 
step-size. 

We use a radial trapping frequency uj r = 2ir x 250Hz, 
axial trapping frequency of u> z = 2n x 50Hz, roughly 
5 x 10 5 atoms, and set the temperature to 1/iK. The 
initial distribution function in position and momentum 

is given by the equilibrium distribution e~^(fc + 2 mw J z ). 



B. Ferromagnetic gas in the x-direction 

The spin density matrix for a ferromagnet pointing in 
the x— direction is: F = ip*ip, where ip = {\: ■ 
Numerically integrating ([5]), we find that for a magnetic 
field of 0.1G (q ~ 15Hz), in the absence of collisional 
relaxation, the populations of the three hypcrfine sub- 
levels coherently oscillate about their equilibrium values, 
as shown in Fig. [5ja). The amplitude of the oscillations 
depends on the magnitude of the spin dependent contact 
interaction. The evolution is not simply a precession, 
rather it involves an oscillation between ferromagnetism 
and ncmaticity. 




FIG. 5: (a) Evolution of central density of itif — component 
(normalized to initial density no = n oo + n u + n~i—i(t = 
0)) as a function of time for negative C2 ( 87 Rb parameters 
(thick line)) and positive C2 ( 23 Na parameters (dashed curve)) 
for q ~ 15Hz. (B ~ 0.1G). The line at n 00 /n = 0.5 is 
a guide to the eye. (b) Thick curve: evolution of wif =0 
population in 23 Na for larger q ~ 60Hz (B ~ 0.3G). Dashed 
curve: suppression of oscillations due to finite relaxation rate 
(r e ; ~ 80ms) in 23 Na within a relaxation time approximation. 
Parameters used in the simulations: ui r — 2n x 250Hz, ui z = 
2tt x 50Hz and N = 5 x 10 5 . 



For C2 < 0(> 0) the population of the mp = sublevel 
oscillates between 0.5 and a larger (smaller) value. The 
frequency and amplitude of the oscillations also depend 
strongly on the magnitude of the quadratic Zeeman en- 
ergy. Increasing q increases the frequency of oscillations, 
but decreases the amplitude (Fig. [SJb)Y. Similar results 
have been obtained by Chang et al. [28| in condensed 
87 Rb, although they considered a different initial spin 
configuration. 

The coherences oscillate at a frequency set by the ex- 
ternal perturbation, in this case the quadratic Zeeman 
energy. The resulting phase factors can be interpreted 
as rotations in spin space. The spin-dependent contact 
interaction couples spin and density, and oscillations in 
the densities of the three components arc seen. Increas- 
ing g, causes the spin vector to rotate faster, leading to 
a rapid averaging to its initial value over the timescale of 
the spin-dependent contact interaction, reducing the size 
of the effect. 

Collisions:— Experiments will also involve collisions 
between the atoms. A detailed derivation of collisional 
relaxation rates can be found in (l6j . We are concerned 
here with the experimental feasibility, hence we use a 
simple relaxation time approximation. The relaxation 



rate is proportional ( 



ao+2a 2 \2 



The amplitude of 



the oscillations in the populations of the three sub-levels 
is directly proportional to the size of c^. 

From the trap dimensions, particle number, and the 
scattering lengths for 87 Rb (|c2 |/co ~ 0.005), we esti- 
mate T e i ~ 20ms. On this timescale, virtually no oscil- 
lations are seen. The story is different for 23 Na where 
the smaller scattering length implies relaxation times al- 
most 4 times longer than in Rb, and the ratio of | C2 1 /co 
is 10 times larger. As shown in Fig.[5]Ja,b) an experiment 
with sodium atoms should be able to detect population 
dynamics. 
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C. Polar state instability 

(b) 




(c) 
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FIG. 6: (a) Collisionless evolution of the tuf = (+1) 
densities (solid/dotted curves) for C2 = — Co, for an initial 
state {^TfA"^} (spin fluctuation), where e = 0.01. The 
dashed curve shows the mp = populations when ijumtiai = 
{-^=,1,—-^=} (nematic fluctuation). The densities are nor- 
malized to the total initial density no = noo + n-n +n_i_i (t = 
0). (b) A gapped and ungapped (solid blue and green) exist 
for C2 = 0. When C2 — > — Co, (dashed blue and green curves) 
the energy of the gapped mode increases (green dashed) while 
the other mode (blue dashed) becomes unstable. The dashed 
orange curve is the imaginary part of this mode, (c) Collision- 
less evolution of central density of = 0, +1 components 
(dotted/dashed respectively) for negative C2 (02/00 ~ 0.05 
and C2 ~ 2Hz) for q <C \c2\n (B ~ lmG). Conservation of 
S z implies that the tuf = — 1 component evolves in the same 
way as mr = 1. (d) Collisionless evolution of coherences 
- top curve is the no-i component and bottom is the ni_i 
component normalized to the initial total density. The ni-i 
component tends to -0.5, coherence is maintained, and the 
final state is a polar state along x. 

In Sec. 11(B) we showed that the polar state of a gas 
with ferromagnetic interactions has a dynamically unsta- 
ble mode at long wavelengths. Here we explicitly demon- 
strate this instability for a trapped gas with ferromag- 
netic interactions. However we also find that this insta- 
bility is not observable in current experiments. Given 
the scattering lengths of 87 Rb, our simulations indicate 
that the collision times are so short that the coherences 
necessary to drive this instability will never have time to 
develop. 

Additionally we explore details of the instability. We 
verify that the polar state is always unstable for small 
q, but the final state is dependent on the magnitude of 
|c2|/co- For small values of this ratio (~ 10 -3 ), finite 
populations of rrif = ±1 develop, but not all the atoms 
are transferred out of the mp = state. Increasing c-i 
decreases the timescale for the onset of this instability, 
as well as increasing the fraction of the population trans- 
ferred. The mp = state is stabilized at large Zeeman 
fields. Note that S z is always a conserved quantity. 

We start with a polar gas with small coherences, which 



seed the instability. We consider seeds with two differ- 
ent symmetries: F = ip*tp where ip = {ipx, tpo, 1} = 
{72' ~75 } ( ncmatic ) or {^f ) !> 75} ( s P in )> where e < 
1, and we only keep terms 0(e) in the density matrix. 
The Wigner functions are assumed to have a Maxwellian 
form in phase space. 

First consider q = dynamics. As the polar state is 
unstable to spin fluctuations, an initial spin fluctuation 
grows exponentially, while a nematic fluctuation merely 
changes the size of the nematic director, but does not 
make it unstable (Fig. E£a)). To enhance the size of the 
effect, we pick C2 = — Co, but we observe an exponentially 
growing spin mode for typical values of c-i as well. 

Although the linear stability analysis predicts an in- 
stability, it does not determine the final state, which de- 
pends in a complicated way on the long time dynamics. 
In Fig. [Ha) we show the collisionless evolution of the 
mp = 0(+l) populations (solid/dashed) starting in the 
state {^j! 1, with e = 0.01. The dashed curve shows 
that for an initial state with nematic order {-^, 1, — 75 }; 
no dynamics is seen. 

In Fig. (6^b) we plot the modes of the polar gas after 
solving |8]). As before there are two modes at C2 = 0. 
One mode has a gap of (c — c 2 )n (green curves), which 
increases as C2 — > — cq. The other ungapped mode (blue 
curves) develops an imaginary part (orange) which in- 
creases as C2 is made more negative. 

The quadratic Zeeman effect stabilizes the polar state. 
Given an initial spin fluctuation about the unmagnetized 
state, this spin precesses in the Zeeman field at a rate 
proportional to q. If the initial state only has a nematic 
perturbation, then a finite q is needed to produce the 
spin fluctuations required to drive an instability. Upon 
increasing q, there arc two effects. First there appear 
wave- vectors where the frequency is real ((5J) , and oscilla- 
tions are observed along with an exponential growth of 
spin populations. Second, the spin vector precesses more 
rapidly, and averages to zero on a timescale of ~ \c2\n, 
ultimately stabilizing the polar state. 

In Fig. |Bfc) we plot the evolution of the central den- 
sity of the mp = 0, +1 sub-levels in time for q ~ lOmHz, 
in a collisionless gas. We take C2n ~ 2Hz, such that 
\c 2 \/c - 0.05, instead of the 0.005 in 87 Rb. In Fig.^d) 
the coherences are shown. The timescale for the onset of 
the instability, as well as the final state depends on the 
size of the off-diagonal seed, and the nature of the pertur- 
bation i.e. whether it is a spin or a nematic fluctuation or 
both, and the strength of the interactions. Here we have 
shown one final state, ipf = {77J, 0, — 775}- Note that the 
final state is not simply an incoherent mixture, unlike 
in a collisional gas. The distribution functions for the 
mp = ±1 components are identical, and no local magne- 
tization is observed. Local structures may be observed 
in an experiment where the trapping potentials are made 
spin dependent. Performing the simulation with a relax- 
ationn rate approximation, taking into acccount collision 
times (r ~ 50ms), and scattering lengths for 87 Rb, we do 
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not observe any dynamics in the polar state. 




1.5 

t(s) 

FIG. 7: Evolution of central density of tuf = component 
(normalized to initial density no = noo + nil + n_i_i(f = 
0)) as a function of time for C2 = 0.5co without (solid) and 
with the relaxation rate approximation (dashed). The initial 
state is F = tp*^p, where -0 = {-^g, 1, -^5} (spin fluctuation), 
where e = 0.01. Trap parameters: 5 x 10 5 atoms, ui r = 2ir x 
250Hz, uj, — 2n x 50Hz. The relaxation time r ~ 50ms. 
For definiteness we assumed 87 Rb atoms, even though the 
scattering lengths in Rb, do not obey this relationship. Owing 
to the strong interactions, one should be able to observe this 
instability even for reasonable collision rates. 

Large anti-ferromagnetic interactions: Next we con- 
sider the case where the spin-dependent and spin- 
independent contact interactions are comparable in mag- 
nitude, and both positive. In Fig. [7] we plot the central 
densities of the mp = atoms as a function of time for 
C2 = 0.5co without (solid curve) and with (dashed) the 
relaxation time approximation, demonstrating the insta- 
bility. Due to the large interaction energies, the dynamics 
is complicated. 

As Endo and Nikuni fiH | have shown, one must be 



careful while considering the relaxation time approxima- 
tion in the c ~ C2 limit where more complicated spin 
dependent collisions may become important. In order to 
fully model the physics, more than one "relaxation-rate" 
may be required. 

V. SUMMARY AND OUTLOOK 

We have addressed the role of the spin-dependent in- 
teraction in the collisionless dynamics of a thermal spinor 
gas. By calculating spin and quadrupolar modes, about 
the magnetized and unmagnetized states, we have shown 
that the normal state of a spinor has a rich array of spin 
excitations compared to its well-studied pseudo-spin i 
counterpart. We numerically calculated the dynamics 
of a non-stationary initial state, finding that the spin- 
dependent contact interaction drives population dynam- 
ics. Finally we provide an explicit demonstration of the 
instability of the polar state to spin fluctuations. 

We conclude that many interesting experiments may 
be done in the normal state of a spinor gas. In order to 
observe some of the physics described here, it will be nec- 
essary to attain a limit where the spin dependent and in- 
dependent interactions are commensurate in magnitude. 
We hope that this work motivates experiments in that 
direction. 
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